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l Under  this  program,  we  have  successfully  developed  the  ultrasonic  impact 
j grinding  technique  for  the  fabrication  of  aspherlcal  geodesic  lens  depressions 
; The  lens  radius  Is  4 mn  and  the  depth  Is  1.374  mm.  15%  of  the  lens  aperture 
I Is  allocated  for  a curved  rounding  at  the  edge  of  the  depression.  The 
’ designed  focal  length  Is  10.2  mm  for  an  aberration  free  aperture  of  6.8  mm. 

| The  technique  of  ultrasonic  Impact  grinding  Is  proved  to  be  precise,  fast 
and  economic.  Reproducible  results  can  be  easily  obtained.  The  tool  wear 
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20.  ABSTRACT  (cont'd.) 

problem  can  be  neglected  If  more  than  two  tools  are  used  for  grinding 
the  lens.  We  have  demonstrated  that  depth  variation  within  7.5/jipfor 
5 depressions  can  be  obtained  using  two  tools.  One  of  the  5 depressions 
has  the  exact  designed  depth.  The  consistency  of  the  lens  depressions 
can  be  further  Improved  If  a third  tool  Is  used.  The  grinding  process 
Is  fast  using  the  200  W Bullen  machine.  It  take^,  only  45  sec.  to  grind 
a 1.4  mm  deep,  8 mm  diameter  depression  In  L1Nb03  substrate.  

Aberration  free  geodesic  lenses  have  been  fabricated  In  T1:L1Nb0, 
waveguides.  The  measured  focal  length  Is  fairly  close  to  the  J 
designed  10.2  mm  even  though  the  lens  profile  Is  slightly  off  the  Ideal 
profile.  The  major  fabrication  problem  encountered  Is  the  polishing  of 
the  'lenses.  There  are  many  fine  scratches  at  the  edge  of  the  depression, 
causing  excessive  propagation  loss  and  In-plane  scattering.  The  measured 
3 dB  focal  spot  size  Is  1.6  pm  as  compared  to  the  theoretical  value  of 
0.32  pm.  The  discrepancy  Is  believed  to  be  the  sum  of  the  following  three 
effects.  (1)  The  resolution  of  a perfect  lOx  microscope  objective  with 
NA  * 0.3  Is  0.93  pm.  (2)  The  use  of  prism  coupler  to  extract  guided  light 
out  of  waveguide  causes  aberration.  The  rutile  prism  Is  uniaxial  and 
Introduces  additional  aberration.  (3)  The  actual  lens  profile  Is  slightly 
off  the  Ideal  designed  profile. 

A ray  tracing  through  geodesic  lens  depression  Is  developed  from  which  the 
geometric  spot  size  and  focal  length  are  calculated.  For  an  f/1.5  lens 
In  L1Nb03,  the  full  lobe  width  of  a Gaussian  beam  truncated  at  the  1/e 
points  Is  about  1.4  pm  assuming  the  free  space  wavelength  is  0.8  pm.  For 
a 1 cm  diameter  lens,  a depth  change  of  IX  can  still  maintain  Its 
diffraction  limited  performance  If  90X  of  the  original  useful  aperture 
Is  utilized.  On  the  other  hand,  the  focal  length  Is  very  sensitive  to 
the  depth  variation.  We  have  shown  that  for  each  unit  decrease  In  total 
depth  for  tool  wear,  there  Is  approximately  a 13  unit  change  In  focal 
length.  With  surface  removal,  the  change  In  focal  length  Is  about  3 units 
per  unit  change  In  depth.  Unless  one  can  control  the  depth  to  better 
than  1.0  pm,  the  Integration  of  geodesic  lenses  In  Integrated  optic 
circuits  Is  an  extremely  difficult  job. 
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INTRODUCTION 


SECTION  I 


Coherent  optical  processing  has  been  considered  for  some  time  to  be 
an  Ideal  technique  for  several  signal -processing  applications  (e.g., 
correlations,  spectral  analysis).  This  Is  due  to  its  unique  capabilities 
of  parallel  processing  and  speed-of-light  operation.  However,  the  size 
of  the  device  and  its  environmental  vulnerability  (to  vibration, 
temperature  variations,  etc.)  precluded  the  field  use  of  the  method  and 
it  remained  a laboratory  tool. 

The  advent  of  integrated  optical  circuits  has  generated  renewed 

interest  in  optical  signal  processing  because  of  their  ruggedness,  small 

size,  and  low  fabrication  cost.  Most  previous  objections  to  optical  processors 

can  be  removed  bv  Integrated  optical  circuits.  In  addition,  for  systems 

requiring  fast  processing  large  data  bandwidths  (larger  than  1 GHzl,  integrated 

optical  circuits  are  prime  contenders.  An  example  of  such  a circuit  is  the 

integrated  optic  spectrum  analyzer,  which  utilizes  the  Bragg  interaction 

between  a surface-guided  coherent  optical  wave  and  a surface  acoustic  wave 

(SAW).  * The  operating  principle  of  the  circuit  is  to  spatially  modulate 

a guided  beam  with  a frequency-modulated  SAW,  then  to  use  a Fourier  transform 

lens  to  focus  the  deflected  beam  on  to  a detector  array.  The  focal 

position  of  the  deflected  beam  is  determined  by  the  angular  deflection 

(or  SAW  frequency).  To  obtain  high  frequency  resolution  from  these 

operations  requires  the  ability  to  perform  the  two  basic  optical  operations, 

collimation  and  focusing,  on  guided  beams  in  optical  waveguides.  Both  are 

limited  by  the  ability  to  construct  aberration-free  single  mode  waveguide 

lenses.  Several  recent  published  papers  have  dealt  with  various  aspects 

3-8 

of  beam  focusing  in  optical  waveguides.  Tn  addition,  methods  of 

3 

aberration  correction  have  been  proposed. 
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This  final  report  summarizes  the  work  performed  during  the  period 
25  April  1978  to  24  April  1979  on  the  contract  ( N001 73-78-C-01 02)  entitled 
"Geodesic  Lenses  for  Integrated  Optic  Circuits".  The  key  goals  of  this 
program  are  (a)  to  develop  and  fabricate  high  precision  aspherlcal  geodesic 
lenses.  The  lens  should  have  diffraction-limited  performance;  (b)  to 
perform  a simulated  computer  error  analysis.  The  effect  of  deviation 
of  the  actual  profile  from  the  Ideal  profile  on  the  lens  performance  will 
be  investigated. 

Under  this  program,  the  ultrasonic  impact  grinding  technique  was 
successfully  developed  to  fabricate  In  Ti  diffused  LINbOg  aberration  corrected, 
single  mode  f/1.5  waveguide  lenses  with  a focal  lens  of  10  mm.  The  program 
goal  is  to  control  the  depth  of  the  geodesic  to  within  1 urn  or  better. 

Other  fabrication  techniques  Investigated  Include  diamond  turning  and 
diamond  generating.  An  extensive  computer  error  analysis  was  performed 
to  determine  the  effects  of  deviation  of  the  actual  machined  profile  from 
the  Ideal  profile  on  the  spot  size  and  focal  length.  A ray  tracing  through 
geodesic  depression  was  developed  from  which  the  geometric  spot  size  and 
focal  length  are  calculated. 

This  report  consists  of  5 sections.  Section  II  discusses  the  design 
of  lens  profile  for  an  aberration-free  focusing  lens.  The  profile  for  an 
f/1.5,  10.2  n*n  focal  length  geodesic  lens  was  computed  for  experiment. 

Section  III  describes  the  development  of  lens  fabrication  techniques,  ultrasonic 
Impact  grinding,  diamond  turning  and  diamond  generating.  Section  IV 
discusses  the  fabrication  and  evaluation  of  geodesic  lenses  fabricated  In 
Ti  diffused  LINbOg.  Section  V reports  the  computer  error  analysis.  In 
particular,  the  utilization  of  ray  tracing  technique  to  determine  the  spot 
size  and  aberration.  Section  VI  reports  conclusions  and  suggestions  for 
future  studies. 
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SECTION  II 


LENS  DESIGN 

Waveguide  geodesic  lenses  are  waveguide-covered  depressions  or 
protrusions.  As  Illustrated  In  Fig.  1,  the  curvature  generates  a 
difference  In  geometric  path  length  for  rays  traversing  different 
portions  of  the  structure;  this  produces  a focusing  effect.  The  optical 
path  length  at  the  center  of  the  depression  Is  longer  than  at  the  sides. 
Because  the  Index  of  refraction  Is  the  same  at  all  points,  the  phase 
fronts  are  curved  toward  the  lens  axis  beyond  the  depression.  In  other 
words,  waveguide  geodesic  lenses  are  always  positive  lenses.  If  the 
depression  Is  of  a particular  aspheric  shape  or  If  the  wavefront  emerging 
from  the  line  Is  appropriately  compensated  by  external  overlayers,  the 
focusing  will  be  aberration  free  (except  for  field  curvature).  The 
aspheric  lens,  the  topic  of  this  report,  Is  formed  by  grinding  and 
polishing  a surface  to  the  aspheric  shape  that  will  give  aberration-free 
focusing  for  the  wavegulded  light.  In  geodesic  lenses,  a rounded  edge 
Is  Included  to  prevent  undesirable  scattering  between  the  planar  guide 
and  the  guide  in  the  lens  depression. 

(a)  General  Design  Principle 

We  use  a geodesic  lens  design  procedure  described  by  Kassal  and 
q 

Marom;  however,  their  procedure  for  designing  the  rounding  has  been 
modified.  This  section  summarizes  the  modified  design  procedure,  which 
Is  based  on  the  principle  of  equivalent  lenses,  and  presents  a generalized 
method  for  calculating  a rounding  with  a smooth  transition  to  the  geodesic. 
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SHAPED  DEPRESSION 


The  design  procedure  entails  the  conversion  of  a rounding  function  to 
a refractive  Index  profile.  The  refractive  Index  distribution  In  the 
remainder  (or  the  Inner  portion)  of  the  lens  Is  calculated  to  achieve 
aberration- free  focusing.  This  refractive  Index  lens  Is  a generalized 
Luneberg  lens.  Finally,  the  refractive  Index  profile  Is  transformed 
Into  an  aspheric  geodesic  profile.  The  physical  properties  of  these 
two  lenses  are  shown  In  Fig.  2. 

Two  lenses  are  equivalent  If  their  Indices  and  surface  profiles 
satisfy  the  condition  of  equal  optical  path  length, where  the  optical  or 
ray  path  Is  determined  by  Fermat's  principle.  Using  this  equivalence. 

It  can  be  shown  that  the  equations  for  transforming  a circularly  symmetric 
flat  guide  with  a refractive  index  distribution  and  a uniform  refractive 
Index  geodesic  lens  profile  are 

n6s  = n6? 

np  = np  (1) 

6 0 = 6F 

where 

n,  rT  are  the  refractive  Indices 

p,  p”  are  the  distances  from  the  axis  of  symmetry 

9,  F are  the  angular  coordinates  with  respect  to  the  optical  axis 

s,  s’  are  the  arc  lengths  of  the  structure  along  a meridional 
cross  section. 

For  planar  structures,  F equals  s.  Throughout  thl s report , barred 
coordinates  and  variables  are  used  for  the  generalized  Luneberg  profile 
and  unbarred  ones  for  the  equivalent  geodesic  lens. 


7034  2*1 


GEODESIC  LENS 


LUNEBURG  LENS 
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It  Is  necessary  to  Include  the  rounded  edge  In  the  design.  The 
rounded  portion  of  the  lens  Is  not  expected  to  have  the  focusing  properties 
of  the  remainder  of  the  lens.  However,  the  effect  of  the  rounded  edge  Is 
taken  Into  account  in  the  lens  profile  calculation.  The  end  result  Is  a 
lens  that  (with  the  rounding  present)  Is  aberration  free  with  the  exception 
of  field  curvature.  The  lens  profile  described  In  Ref.  9 (in  which 
a circular  rounding  was  used) were  not  sufficiently  smooth  to  prevent 
excessive  scattering  at  the  lens  rounding  transitions.  This  was  due  to  a 
rapid  curvature  change,  which  was  necessary  to  match  the  slope  at  the 
transition.  This  rapid  curvature  change  was  also  accompanied  by  a dis- 
continuous second  derivative  in  the  profile  functional  shape.  A generalized 
rounding  design  procedure  is  described  below. 

The  rounding  is  shown  schematically  In  Fig.  3.  Several  constraints 
are  placed  on  the  rounding  function,  r,  to  reduce  the  scattering  and  the 

.] 

insertion  loss  associated  with  the  disturbance  Introduced  by  the  lens  in 
the  waveguide: 

(I)  The  slope  of  r at  the  lens  edge  Is  zero 

(II)  The  rounding  function  r and  the  transition  between  the 
rounding  and  the  geodesic  profile  are  continuous 

(III)  The  slope  of  r and  the  slope  of  the  geodesic  profile  are 
continuous  at  the  transition 

( 1 v ) The  second  derivative  of  r and  the  second  derivative  of 
geodesic  profile  are  continuous  at  the  transition. 

To  design  a Luneberg  lens,  the  proposed  rounding  function,  r(x),  is 
transformed,  using  Eq.  1,  Into  a refractive  Index  profile.  The  algorithm 
for  performing  the  transformation  Is  carried  out  digitally  by  stepping 
the  variable  x = R-p  from  0 to  the  rounding  edge  (l.e.,  x.+^  = x^  + Ax), 
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where  R is  the  lens  radius  and  the  subscript  i refers  to  the  annulus 
number.  For  each  Ax,  the  corresponding  As  is  given  by 


Rewriting  Eq.  1 in  different  form  yields 

nAs  = n(xi+1)  Ax  = n(x1+1)(x1+1  - x^)  (3) 

n(xi+1)  - np/p  = n(R  - xi+1)/(R  - x1+]),  (4) 

where  - x^  is  the  incremental  Ax  corresponding  to  Ax,  x^  is  the  previous 
value  of  x obtained  in  the  preceding  annulus,  and  x^  and  W(x..+j)  are  the  values 
to  be  calculated.  Eliminating  n(x^+-|)  from  Eqs.  3 and  4 gives 

RAs  + Xj  (R  - x,+1) 

7 = ] 1±!_  (5) 

xi+l  R - x i+1  + As  • 

Substituting  this  result  Into  Eq.  4 determines  ri(x^+1). 

To  determine  the  function  r that  satisfies  conditions  (iii)  and 
(iv)  above  requires  the  tracing  of  rays  through  the  system.  The  initial 
problem  associated  with  ray  tracing  is  to  find  which  ray  coming  from  the  focal 
point  will  penetrate  a specified  distance  (i.e. , x = x^)  into  the  lens. 

Light  rays  are  known  to  satisfy  the  following  condition  on  a surface  of 
revol ution: 


p sin  a * h,  (6) 

where  a is  the  angle  between  the  light  path  and  the  meridian  at  p,  and  h is 
a constant  equal  to  the  shortest  distance  between  the  light  trajectory  on 
the  geodesic  surface  and  the  axis  of  rotation.  Thus,  for  any  x,  the 
angles  a and  (Figure  4)  are  given  by 


(7) 

(8) 
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Since  the  mapping  between  the  generalized  Luneberg  lens  and  the 

g 

uniform  refractive  Index  geodesic  lens  Is  conformal,  the  angles 
a and  ty  are  the  angles  required  for  ray  Initialization. 

After  calculating  the  refractive  Index  distribution  of  the  Luneberg 
lens  and  transforming  It  to  create  a constant-refractive-index  geodesic 
lens,  there  Is  still  the  problem  of  satisfying  conditions  (111)  and  ( 1 v) . 

This  Is  accomplished  by  varying  the  rounding  function f . For  example.  In 
the  function 

r = Ax2  + Bx3,  (9) 

the  coefficients  A and  B are  varied.  The  points  corresponding  to  x - xt  are 
calculated  based  on  the  focusing  requirement;  the  points  x < xfc  are  calculated 
using  Eqs.  2 through  5.  The  slope  and  second  derivative  are  calculated  using 
a two-  and  three-point  method,  respectively,  at  the  transition  between  the 
rounding  and  the  lens  region.  The  slope  calculated  using  the  points  x = - Ax 

and  x = xt  is  required  to  equal  r '(xt  - Ax/2).  The  second  derivative  using  the 
points  x = xt  - Ax,  xt,  xt  + Ax  Is  required  to  equal  r"(xt).  After  finding 
values  of  A and  B that  satisfy  conditions  (111)  and  (1v),  the  design  of  the 
whole  lens  can  be  completed  by  Incrementing  x ort  and  requiring  that  rays 
originating  at  the  desired  focal  point  exit  parallel  to  the  optical  axis. 

Each  ray  at  an  angle  - Ai|>  * or  x^+^  = x^  + Ax  determines  the  Index 
of  the  central-most  annulus  through  which  It  transits. 

Two  important  properties  of  the  geodesic  lens  have  not  been  discussed 
In  the  literature.  First,  the  geodesic  lens  has  no  chromatic  aberrations 
because  of  the  uniform  refractive  Index.  This  Is  Independent  of  the  fact 
that  the  material  refractive  index  Is  a function  of  the  wavelength.  Second, 
the  lens  Is  mode  independent.  This  may  allow  multi-  and  single-mode 
optical  systems. 


(b)  Design  of  f/1/5  lens 


The  profile  of  the  surface  depression  for  an  aberration-free 
geodesic  lens  was  derived  using  the  design  principles  outlined  above. 

An  analysis  of  radiation  In  curved  waveguides  has  shown  that  It  Is 
necessary  to  Include  In  the  design  of  the  rounding  a constraint  on  the 
curvature  of  the  lens  profile.10  In  the  design  of  an  f/1.5  lens  for 
this  contract,  it  was  necessary  to  use  an  aberration-free  aperture  of 
only  85%  of  the  total  lens  diameter  to  avoid  cut-off  at  every  point 
on  the  geodesic.  In  view  of  this  experience  we  have  Included  In  our 
design  program  "DILIGENT"  an  algorithm  to  check  that  the  radius  of 
curvature  of  the  geodesic  always  exceeds  a minimum  curvature  below  which 
the  waveguide  is  beyond  cutoff.  Figure  5 shows  the  profile  computed 
for  an  f/1.5,  10.2  mm  focal  length,  geodesic  lens.  The  lens  diameter  is 
8 mm,  however,  the  useful  (aberration-free)  aperture  is  6.8  tim.  The 
depth  at  the  depression  center  is  about  1.37  mm. 

(c)  Radiation  at  the  bottom  of  lens  depression 

We  have  studied  the  radiation  from  the  bottom  of  the  geodesic  lens 
into  the  substrate.  The  derivation  of  fiber  leaky  waves  by  Snyder  and 
Love11  was  used  here.  Figure  6 shows  a guided  wave  incident  onto  a curved 
waveguide-substrate  interface.  In  order  to  match  the  field  at  the  boundary 
part  of  the  energy  will  be  radiated  into  the  substrate.  The  field  in  the 
substrate  first  decays  exponentially  and  then  starts  to  oscillate  at  a 
point  xt,  determined  by 


*t  ■ f - 


sin2e 


sin  6 
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5.  Geodesic  lens  profile  for  f/1.5  lens 
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where  p Is  the  radius  of  curvature  of  the  Interface,  0 Is  the  critical 
angle  between  the  nj  and  n^  media,  and  0 1s  the  corresponding  angle  of  the 
guided  mode.  The  decay  constant  a for  the  field  In  the  substrate  Is 
approximately 

a «*4.2052  - 2.202  ^ 

- 0.148  ^ 

for  a simplified  model  of  step  Index  waveguide.  For  a Tl-dlffused  LiNbOg 
waveguide,  we  obtain 


xt  ■ | x 0.0043  for  0C  = 84,55°  and  0 = 86.05°  (neff  ^2.205 

Excessive  radiation  occurs  If  xt  Is  smaller  than  the  decay  length  of  the 
evanescent  field  (1/a)  In  the  substrate.  Thus  radiation  loss  is  negligible 
If 


which  means 


x*  >- 
t a 


p >0.44  mm 


An  alternate  derivation,  based  on  a different  approximation,  finds 


‘t  * f (1  - 


which  finally  leads  to  a minimal  radius  of  curvature  of  0.88  mm,  much 
less  than  the  radius  designed  for  most  practical  lenses.  We  conclude 
that  radiation  from  the  bottom  of  our  geodesic  lenses  into  the  substrate 
Is  not  substantial  and  can  be  neglected. 


SECTION  III 


FABRICATION  OF  GEODESIC  LENS  DEPRESSIONS 
(a)  Ultrasonic  Impact  Grinding 

Ultrasonic  Impact  grinding  Is  a complicated  process  governed  by 
many  parameters;  namely,  tool  shape,  tool  material,  grinding  pressure, 
stroke  amplitude,  workpiece  material,  slurry  speed,  abrasive  concentra- 
tion and  grit  size.  The  operational  conditions  under  which  the  grinding 
experiment  was  performed  are  briefly  described  below. 

(1)  Tool  shape 

There  are  many  different  shapes  of  the  tool  holder  which  can 
be  used  In  the  ultrasonic  Impact  grinding  machine.  Tapered  tools  with 
the  bottom  cross  section  smaller  than  the  top  have  the  advantage  of 
amplitude  amplification.  As  a result,  the  cutting  speed  Is  increased. 
However,  tapered  tools  do  not  have  long  term  stability,  a constant 
tuning  of  the  resonance  condition  Is  necessary  for  precision  work. 

On  the  other  hand,  straight  cylindrical  tools  have  slower  cutting  rate 
but  stay  In  resonance  for  a long  time.  The  designed  tool  holder  In  the 
experiment  has  a straight  shape.  It  takes  an  average  45  minutes  to  grind 
an  8 mm  diameter,  1.5  mm  deep  depression  In  LINbO^  using  20  pm  Al^O^ 
abrasive  slurry. 

(2)  Tool  material 

Tool  material  should  be  selected  based  on  Its  Impact  resistance 
to  wear.  Good  tool  wear  is  obtained  with  material  which  has  high  strength 
and  good  ductility.  The  tool  tips  we  used  are  made  of  cold  rolled  stainless 
steel  420.  The  material  was  annealed  before  machining. 

(3)  Abrasives 

Abrasive  grains  are  the  cutting  edge  of  the  tool.  They  are 
dispersed,  with  the  help  of  suspension  agents.  In  the  water  that  flows 
between  the  workpiece  and  the  tool.  Conmonly  used  abrasives  are  boron 
carbide,  silicon  carbide,  and  aluminum  oxide.  We  choose  20  pm  AlgO-j  for 
rough  grinding  and  9 pm  A1203  for  fine  grinding.  The  abrasives  are  mixed 
with  water  at  a ratio  of  about  1:1  In  volume.  We  have  also  purchased  an 
additional  recirculation  pump  for  3 pm  A1203  powder.  The  pump  Is 
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used  to  study  the  Improvement  of  surface  roughness  when  finer  grits 
are  used. 

During  this  program,  we  have  designed  Impact  grinding  tools  for 
f/1.5  geodesic  lenses.  The  lens  radius  Is  4 nm  and  the  depth  Is  1.374  mm. 

15%  of  the  lens  aperture  Is  allocated  for  a curved  rounding  at  the  edge 
of  the  depression.  The  designed  focal  length  Is  10.2  mm  for  an  aberration 
free  aperture  of  6.8  mm.  A survey  of  local  machine  shops  revealed  that  a 
5 ym  tolerance  In  tool  tip,  which  bears  the  designed  lens  profile.  Is 
about  the  limit  that  an  experienced  machinist  can  achieve  with  conventional 
machine  shop  equipment.  To  reduce  the  tolerance  to  less  than  5 ym  requires 
different  fabrication  techniques  which  are  not  readily  available  In  most  of 
the  machine  shops.*  We  have  purchased  the  tool  tips  with  5 ym  tolerance 
for  the  f/1.5,  F * 10.2  urn  lens  experiment.  The  blueprint  of  the  tool  tip 
design  Is  shown  In  Fig.  7.  The  tools  are  designed  for  our  Raytheon  100  W 
ultrasonic  Impact  grinding  machine.  However,  the  tool  tip  Is  also  adaptable 
to  the  new  high  precision  Bullen  machine. 

It  has  been  reported  that  ultrasonic  Impact  grinding  technique  Is 
capable  of  fabricating  depressions  complementary  to  the  tool  shape  with 

I 

1/4  ym  tolerance.  The  operation  requires  more  than  one  tool  to  avoid  the 
tool  wear  effect.  It  Is  estimated  that  three  different  tools  are  needed 
to  obtain  precise  lens  depressions  with  tolerance  of  1 ym  or  better.  Using 
the  tools  purchased  from  the  local  machine  shop,  we  have  ground  11  geodesic 
lenses  In  LINbO^  substrates.  Two  tool  tlpswereused  to  fabricate  each  lens 
depression,  one  for  rough  grinding  and  one  for  fine  grinding.  The  results 
of  grinding  these  lenses  are  now  presented. 

*We  are  unable  to  find  a local  machine  shop  which  Is  capable  of  fabricating 
tools  with  1 ym  profile  tolerance.  On  the  other  hand,  there  are  techniques 
in  the  field  of  optical  component  fabrication  where  the  dimensional 
tolerances  on  the  order  of  optical  wavelengths  are  frequently  achieved. 

The  optical  techniques  are  diamond  turning  and  diamond  generating.  The 
diamond  turning  machine  Is  very  expensive  and  the  turning  technique  does 
not  work  for  stainless  steel.  The  diamond  generating  seems  a viable 
technique  and  Is  less  expensive. 
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Figure  7.  Design  of  the  tool  tip  for  ultrasonic  impact  grinding 


(b)  Experimental  Results  of  Ultrasonic  Impact  Grinding 

Me  have  fabricated  11  geodesic  lens  depressions  In  LiNbOg 
substrates  using  different  combinations  of  tool  tips.  The  order  of  using 
tools  and  the  depression  depth  are  summarized  in  Table  1.  Due  to  the  lack 
of  depth  and  position  control  from  the  Raytheon  machine  all  the  lenses  are 
ground  deeper  Into  the  substrate  than  the  desired  depth.  The  overall 
depth  of  each  depression  measured  from  the  top  surface  (true  lens  depth  + 
over  grind)  is  also  listed  In  the  Table. 

Tool  tip  No.  2 was  used  to  rough  grind  LiNbO^  samples  2A,  2B  and  2C. 

Sample  2A  had  an  oval  shape  because  the  sample  moved  during  the  rough 
grinding  process.  After  this  run,  soft  wax  was  used  all  the  time  to  keep 
the  sample  in  position.  When  the  sample  2B  was  rough  ground,  the  tool 
alignment  mechanism  was  out  of  order,  the  resulting  depression  was  tilted 
at  angle  more  than  10  degrees. 

Sample  2C  was  broken  when  we  tried  to  tune  the  tool  into  resonance  while  the 
tool  is  in  contact  with  the  sample  surface.  After  losing  three  LiNbOj  samples, 
we  have  gained  enough  experience  to  successfully  grind  the  rest  of  the  8 samples. 
Repositioning  the  sample  for  fine  grinding  requires  some  skill,  particularly 

with  the  Raytheon  machine  which  is  not  equipped  with  precision  translation 
stages.  We  were  very  successful  in  positioning  all  samples  except  sample  4B 
which  was  slightly  misplaced. 

The  conclusions  we  have  reached  so  far  are  discussed  below. 

(1)  The  tool  wear  of  the  stainless  steel  tool  is  slightly  less  than  2%. 

The  degree  of  tool  wear  depends  on  the  parameters  used  in  the 
grinding  process.  More  careful  experiments  are  required  to  establish 
the  optimum  operational  condition  for  a minimum  tool  wear. 

(2)  The  grinding  rate  is  fast  even  with  100  W Raytheon  machine.  It  takes 
an  average  45  minutes  to  grind  the  depression.  The  cutting  rate  is 
slower  when  the  contacted  working  area  becomes  larger.  For  the  last 
few  thousandths  where  the  tool  tip  is  embedded  in  the  substrate,  the 
cutting  rate  is  about  0.5  mil  per  minute. 

(3)  The  tool  alignment  Is  critical  to  the  fabrication  of  high  precision 
lens.  When  the  tool  is  not  perpendicular  to  the  substrate  surface, 
the  lens  thus  made  will  be  tilted  at  an  angle  equal  to  the  angular 
misalignment.  If  we  specify  that  the  difference  between  two  opposite 
points  on  the  lens  edge  be  smaller  than  1 ym,  the  angular  misalignment 
should  be  less  than  7 x 10  degree  (or  25.8  sec)  assuming  the  lens 


Is  In  the  unit  of  volts.  (1.308  mV  = 1.0  ym)  The  designed  lens  depth  Is  1.797  V. 
Except  sample  3E,  the  depth  variations  of  all  the  samples  after  fine  grind  are 
within  0.3%. 


diameter  Is  8.0  m.  Angular  alignment  of  better  than  10  sec  can  be 
achieved  by  optical  means.  We  have  used  a tool  maker  square  to  check 
the  perpendicularity  before  the  grinding  process  is  started.  The 
angular  accuracy  obtained  is  better  than  0.1°.  To  compensate  for 
the  slight  angular  misalignment,  each  lens  was  ground  deeper  into  the 
substrate  than  the  designed  depth.  Subsequently,  the  top  flat  surface 
was  re-established  using  the  lens  axis  as  the  reference.  The  whole 
substrate  was  tilted  with  respect  to  the  polishing  pad  such  that  the 
polished  flat  surface  touched  the  edge  of  the  depression  simultaneously. 
This  compensation  technique  worked  well  for  substrates  having  only  one 
lens.  We  plan  to  incorporate  an  optical  instrument  into  the  new 
ultrasonic  machine  so  that  the  angular  requirements  are  met  before 
grinding. 

(4)  The  technique  of  ultrasonic  impact  grinding  is  proved  to 
be  precise,  fast  and  economic.  Reproducible  results  can  be  fairly  easily 
obtained.  The  tool  wear  problem  can  be  neglected  if  more  than  two  tools 
are  used  for  grinding  the  lens.  Reference  to  Table  1 reveals  that  the 
five  samples  fine  ground  with  tool  tip  No.  6 have  resulting  lens 
depressions  with  depth  variation  within  10  mV  (^7.5  ym).  The  consistency 
of  these  lens  depressions  will  be  better  than  7.5  ym  if  a third  tool 
tip  is  used. 

The  Raytheon  ultrasonic  impact  grinding  machine  is  not 
designed  for  high  precision  work.  We  have  to  go  through  unnecessary 
steps  to  preserve  the  lens  profile.  Even  with  the  great  difficulties 
encountered,  however,  we  were  able  to  fabricate  depressions  within 
several  ym  of  the  designed  profile.  In  order  to  meet  the  goal  of 
precision  machining,  an  upgraded  Impact  grinding  machine  was  purchased 
using  internal  funds.  The  new  machine  has  the  capability  of  controlling 
the  depth  within  1 ym.  The  preliminary  test  demonstrated  a cutting 
speed  about  50  times  faster  than  the  old  machine.  It  took  only  45  sec 
to  grind  a 1.4  mm  deep,  8 mm  diameter  depression  in  LINbOg  substrate. 

The  surface  finish  was  also  Improved  by  reducing  the  tool  strokes 
during  the  grinding  of  the  last  25  ym. 
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(b)  Diamond  Turning  and  Diamond  Generating 

Task  1 of  the  statement  of  work  of  this  contract  involved  an 
investigation  of  the  use  of  numerically  controlled  precision  machining 
for  formation  of  geodesic  lenses.  Two  techniques,  singie  point  diamond 
turning  and  diamond  generating,  were  investigated. 

Diamond  turning  is  a precision  cutting  process  using  a single 
point  diamond  as  the  cutting  edge  to  form  the  geodesic  lens  depression. 
Diamond  turning  machines  are  commercially  available  from  Moore  Special 
Tool  Co.  for  fabricating  aspherical  optics.  The  major  utility  of  this 
technique  is  for  low  cost,  high  precision  infrared  optics  where  tool 
marks  on  the  surface  can  be  tolerated.  At  these  wavelengths,  the  machined 
surfaces  are  readily  usable  without  additional  polishing.  When  a LiNbOg 
single  crystal  Is  turned,  the  surface  shows  a hazy  appearance  in  areas 
termed  as  unfavorable  crystalline  regions  resulting  from  anisotropic 
crystalline  properties.  A delicate  optical  polish  after  turning  is 
therefore  necessary  for  fabricating  integrated  optic  devices.  This 
observation  agreed  with  our  experimental  results,  5 ym  peak-to-valley  surface 
roughness  was  measured  In  the  unfavorable  regions  of  a flat  diamond  turned 
LiNbOg  substrate. 

Although  the  state-of-the-art  diamond  turning  machines  have  a 1/4 
ym  control  accuracy  In  both  translational  and  rotational  movements,  the 
resulting  lens  profiles  on  LiNbOg  have  depth  variation  of  up  to  4.0  ym 
and  the  profile  tolerance  of  about  3.0  ijji.  The  poor  depth  control  is 
attributed  to  troubles  of  gauging  the  depth  and  can  be  improved.  However, 
the  problem  of  profile  tolerance  of  3.0  ym  is  due  to  the  fact  that  a large 
tool  pressure  Is  required  to  cut  LiNbOg  substrates.  This  is  an  intrinsic 
problem  associated  with  the  anisotropic  characteristics  of  the  material. 

Diamond  generating  employs  a diamond  impregnated  tool  to  grind  the 
surface  depression.  The  grinding  Is  controlled  by  precision  slides  and 
spindles  which  are  Interfaced  with  a computer.  This  technique  has  been 
employed  successfully  to  fabricate  large  diameter  aspherical  lenses 
(3-dimensional  bulky  lens).  Special  tools  need  to  be  developed  for 
small  radius  aspherical  shapes.  A directed  subcontract  for  diamond 
generating  of  4 f/1.5  geodesic  lenses  In  LiNbOg  substrate  has  been  given 
to  a local  optical  shop,  OPTI  Systems,  Inc.  The  profile  tolerances  of 


the  lenses  are  specified  at  1/2  vim  for  two  samples  and  1 urn  for  the  other 
two  samples.  The  experimental  results  from  OPTI  Systems  Indicated  that 
the  surface  damage  due  to  diamond  grinding  was  worse  than  expected, 
pieces  of  LiNbO^  larger  than  the  diamond  grit  size  chipped  off  the 
surface  during  grinding.  The  depth  of  this  subsurface  damage  varied 
from  6 to  10  ym  peak- to-val ley.  We  conclude  that  diamond  generating 
is  not  an  adequate  technique  for  fabricating  aspherical  depressions  in 
LiNbO-j  crystals. 


SECTION  IV 


FABRICATION  AND  EVALUATION  OF  GEODESIC  WAVEGUIDE  LENSES 


Three  LiNbO^  samples,  3A,  4A,  and  4C,  were  polished  during  the 
second  quarter  of  this  program.  The  task  of  maintaining  the  lens  profile 
during  optical  polish  was  a success.  We  have  measured  a depth  change 
less  than  1.0  urn  after  polish  for  all  three  samples.  Optical  waveguides 
were  fabricated  on  these  samples  in  a conventional  way.  Single  mode 
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waveguides  were  formed  by  diffusing  200  A E-beam  evaporated  Ti  at  1000°C. 

The  diffusion  time  was  5 hours.  Although  the  He-Ne  laser  light  could 
propagate  in  the  planar  waveguide  region,  all  the  guided  light  terminated 

O 

right  at  the  lens  edge.  We  had  tried  to  add  200  A more  Ti  onto  samples 
4A  and  3C  and  repeated  the  diffusion  process.  Deeper  waveguides  and  larger 
index  change,  however,  did  not  eliminate  the  radiation  problem.  We 
believe  the  radiation  is  due  to  a sharp  corner  at  the  edge  caused  by  the 
polishing  of  the  top  surface.  Since  there  was  no  depth  control  in  the 
Raytheon  machine,  we  ground  the  depression  all  the  way  into  the  substrate. 

The  overground  portion  was  then  removed  by  removing  the  top  flat  surface. 

Using  the  new  machine  we  should  be  able  to  control  the  grinding  process 
with  high  precision,  thus  eliminating  the  process  of  surface  removal. 

Since  contract  funds  were  not  sufficient  to  grind  and  polish  more 
geodesic  lenses  using  the  high  precision  ultrasonic  impact  grinding  machine 
which  was  installed  and  operational  in  early  February  of  1979  we  decided, 
with  the  approval  from  NRL  (Naval  Research  Laboratory),  to  reprocess 
sample  3A,  4A  and  4C.  The  Ti  diffused  waveguide  layer  was  removed  and 
the  edges  of  the  lens  depressions  were  ground  and  repolished.  The  sharp 
corner  at  the  lens  edge  was  removed  to  ensure  continuous  propagation  of 
waveguided  light.  As  expected,  the  resulting  lens  profile  was  slightly 
deviated  from  the  original  designed  one.  Sample  4A  was  delivered  to  NRL 
at  the  end  of  this  program. 

The  optical  waveguide  was  formed  on  sample  4A  by  diffusing  300A  Ti 
into  the  LiNbO^  substrate  in  an  oxygen  atmosphere.  The  diffusion  temperature 
was  1000°C  and  the  diffusion  tine  was  6 h.  Before  diffusion,  the  Ti  metal 
was  oxidized  at  600°C  for  4 h.  The  waveguide  can  support  two  TE  and  two  TM 
polarization  modes.  However,  only  th«  lowest  order  mode  can  propagate  through 


the  lens  depression.  At  GaAJtAs  laser  wavelength,  the  waveguide  becomes 
single  mode.  For  convenience,  an  He-Ne  laser  was  used  as  the  light  source 
in  the  experiment.  The  focal  length  of  the  geodesic  lens  was  determined 
by  the  deflection  angle  of  the  He-Ne  laser  light  coupled  out  by  a rutile 
prism  positioned  after  the  lens.  The  deflection  angle  was  measured  from 
the  beam  displacement  and  the  distance  between  the  lens  and  the  detector. 

To  reduce  the  effect  of  prism  aberration  on  the  measurement,  the  prism  was 
positioned  such  that  the  focal  point  is  approximated  at  the  hypotenuse  face 
of  the  prism. 

Fig.  8 shows  the  measured  focal  length  normalized  to  the  radius  of 
the  depression  as  a function  of  the  lateral  displacement  of  the  input  beam. 

The  focal  length  remains  fairly  constant  over  85%  of  the  lens  depression 
the  designed  useful  aperture.  The  slight  tendency  for  the  focal  length 
to  become  shorter  as  the  light  beam  moves  away  from  the  center  of  the  lens 
could  be  due  to  either  the  profile  deviation  caused  by  the  additional  polishing 
mentioned  previously,  or  the  experimental  inaccuracy. 

The  intensity  distribution  near  the  focal  plane  was  measured  using  an 
expanded  He-Ne  laser  beam.  The  rutile  prisms  were  used  to  couple  the  laser 
light  into  and  out  of  the  waveguide.  After  the  output  prism,  a relay  lens 
( 1 Ox  microscope  objective  with  NA  = 0.30)  imaged  and  magnified  che  focal 
spot  to  a 7 - pm  slit  in  front  of  the  detector.  New  rutile  prisms  with  60° 
base  angle  were  used  in  the  experiment.  As  a result,  a higher  power  microscope 
objective  with  short  working  distance  can  be  utilized  to  reimage  the  focal  plane. 
An  iris  at  the  input  end  controlled  the  width  of  the  input  beam.  Fig.  9 
shows  the  intensity  distribution  when  the  full  aperture  was  illuminated. 

The  3dB  spot  size  was  1.6  ym  as  compared  to  the  theoretical  spot  size  of 
0.32  ym.  The  discrepancy  is  the  result  of  the  following  three  effects: 

(1)  the  resolution  of  a perfect  lOx  microscope  objective  is  0.93  ym 

(2)  the  presence  of  prism  coupler  introduces  aberration  (3)  the  deviation  of 
lens  profile  from  the  ideal  profile.  Nevertheless,  the  spot  size  is  the 
smallest  ever  measured.  The  side  lobes  to  the  right  of  the  main  peak  could 
be  due  to  severe  in-plane  scattering  or  lens  aberration.  When  the  relay  lens 
was  moved  by  a small  amount  in  the  longitudinal  direction,  the  intensity 
distribution  shown  in  Fig. 10  was  obtained.  At  this  position  the  spot  size  is 
bigger  (3dB  width  2.5  ym),  however,  the  side  lobes  in  Fig. 10  become  diffused. 
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Fig.  8 Measured  focal  distance  vs.  lateral  displacement 
of  Incident  beam  for  geodesic  lens  fabricated 
In  T1:L1Nb0g  waveguide. 


(I 


-26- 


Intensity  distributed  at  the  focal  point  for  the  geodesic 
lens  fabricated  in  a TirLiNbO^  waveguide. 


Fig. 10  Intensity  distribution  near  the  focal  point  for  the 

geodesic  lens.  The  relay  lens  (microscope  objective) 
was  placed  at  a position  where  side  lobe  levels  are  lo 


SECTION  V 


COMPUTER  ERROR  ANALYSIS 


(a)  Ray  Tracing  Theory 

In  this  section  we  summarize  a ray-tracing  technique  for  tracing  rays 
on  a rotatlonally  symmetric  geodesic  lens.  By  Fermat's  principle  It  can 
be  shown  that  light  rays  travel  along  geodesic  lines.  The  differential 
equations  that  the  geodesics  of  a surface  of  revolution  satisfy  are: 


2 d0  d£ 
p ds  ds 
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where  0 Is  the  angle  of  rotation,  p Is  the  distance  from  the  axis  of  rotation, 
s Is  the  path  length  traveled  by  the  light  ray  Inside  the  surface,  and  g(p) 

Is  the  function  that  describes  the  profile  which  generates  the  surface  of 
revolution,  g*  and  g"  are  the  derivatives  of  g(p)  with  respect  to  p.  The 
differential  equations  given  by (10 ) and (11)  can  be  solved  by  separation  of 
variables  and  rewriting  the  equations  in  the  form  of  first-order  linear 
differential  equations.  In  a cylindrical  coordinate  system  the  solutions  of 
the  first  order  differential  equations  are  reduced  to 
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where  0Q  Is  the  initial  azimuth  of  the  entering  light  ray,  pQ  is  the  outer 
radius  of  the  lens,  c = pQ  sin  0Q,  and  g‘  is  the  derivative  of  the  profile 
with  respect  to  the  radius  p of  the  cylindrical  coordinate  system.  e(p)  is 
the  azimuth  and  s(p)  is  the  path  length  of  the  ray  measured  from  the  point  of 
entrance  into  the  geodesic  lens. 
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The  ultimate  goal  of  this  task  Is  to  write  a computer  program  for 
tracing  rays  through  rotatlonally  symmetric  depression  lenses.  In  order 
to  determine  the  propagation  direction  of  the  light  the  tangent  vector 
of  the  path  was  calculated.  In  a rectangular  coordinate  system  the  unit 
tangent  vector,  t,  Is  given  by: 


(P.0) 


1 ♦ (9*)* 


c sin  6-i-t 

- r-]e] 


r-s1n0 


,2-c2 


1 ♦ <9')‘ 


-)'Z  * 


c-31 

L p 


2 2 
P - C 


d + (g 


where  ejs  are  unit  vectors  along  the  vn  axis.  The  path  of  a light  ray 
on  the  rotationally  symnetric  lens  is  symmetric,  such  that  p (the  radius 


used  as  a parameter  to  describe  the  path)  decreases  to  a minimum  of  pQ  sin  80 


and  then  increases  back  to  p„ 


Thus  the  azimuth,  8e»  of  the  point  where  the 


light  ray  departs  from  the  lens  Is: 

sin  8 0 71/9 

0e  ■ eo  * l2c  ! 


P -c 


Substitute  pQ,  0e  for  p,8  in  the  equation  of  the  tangent  vector,  Eq.  (5)  then 
t(p0,8e)def1nes  the  direction  of  the  outcomlng  ray.  This  In  turn  determines 


the  "focal  distance",  f,  of  the  ray.  f Is  defined  as  the  distance  between  the 
lens  center  and  the  point  where  the  outcoming  ray  crosses  the  lens  axis  and 
can  be  expressed  In  the  simple  form 


f I 


pQ  sin  e0 


slnCe0  + V 


The  total  optical  path  length  of  the  ray  Is 


Vin0o 


se  - so  + 

Pr 


1 + fa’)2  1/2 

( 22  > pdpl 

P - c^ 


The  ray- tracing  technique  summarized  above  was  checked  on  results 
which  are  proven  correct.  Two  special  uses  of  the  above  equations  will  be 
outlined  below.  Their  application  to  lens  tolerance  analysis  will  be  pointed 


First  example  Is  the  Investigation  of  the  focusing  properties  of 
spherical  depressions.  For  spherical  depressions 

. _ „ • «o2  Sl»2e„  * C-2  «0)«\  (,8) 


®o  ♦ \j~  arc  sin  (- 


slnV  - R 

o 


and  for  the  limiting  case  of  paraxial  ray 

11m  (—  ) * 2(1  - cosiH 


^ P?) 


where  ^ Is  the  polar  angle  and  R Is  the  radius  of  the  sphere,  a portion  of 
which  creates  the  depression.  This  agrees  with  published  results  for 
spherical  depressions. 

In  tolerance  analysis  of  geodesic  lenses  this  example  of  spherical 
depression  is  useful,  since  the  explicit  and  relatively  simple  equation  for 
the  focal  length  allows  one  to  estimate  changes  In  the  focal  length  of 
geodesic  lenses  resulting  from  slight  deviations  In  the  shape  and  depth  of 
the  lens. 

Another  example  for  rotational ly  symmetric  surfaces  Is  a toroidal  section 
(g(p)  * r - /r2  - (P0  - p)2)  generated  by  the  section  of  a circle  of  radius  r. 
The  Importance  of  this  surface  of  rotation  is  that  it  can  be  viewed  as  a 
rounding  to  the  lens,  and  In  that  sense  It  Is  a prototype  of  any  kind  of 
rounding.  To  know  the  focusing  behavior  of  such  surfaces  is  important  in 
the  construction  of  rounded  edge  geodesic  lenses.  In  fact,  when  a particular 
focal  length  is  desired  In  connection  with  circular  rounding  the  question  is 
whether  any  light  beam  parallel  to  the  central  ray  when  entering  the  toroidal 
section  will  terminate  its  path  In  the  direction  of  the  desired  focal  point. 
The  solution  of  this  problem  Involves  elliptical  Integrals  of  the  first, 

F(X,  q),  and  third,  m(X,  n,  q),  kinds. 


0e  ■ eQ  ♦ |A  F(X,  q)  +B  m(X,  n,  q) 
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The  focal  distance  for  each  ray  which  path  stays  entirely  on  the  toroidal 
rounding  section  is  then  given  by  equation  (7).  These  equations  can  be 
applied  to  determine  the  effects  of  variation  in  the  lens  rounding. 

(b)  Computer  Program 

A ray  trace  program  has  been  developed  using  Eqns.  (14)-(17).  This  was 
not  a trivial  task  in  that  the  available  integration  routines  were  not 
adequately  accurate.  This  was  due  to  the  fact  that  the  integrand  contains 
at  least  one  singularity  and  in  some  cases  two.  In  addition  the  normal 
technique  of  integration  by  parts  to  eliminate  a singularity  of  this  form 
did  not  work  because  the  test  case  profile  has  a discontinuous  g"(p). 

This  test  case  profile  defines  an  aberration-free  lens  with  focal  point  on 
the  lens  edge. 

An  algorithm  for  calculating  the  integral  has  been  found  and  included 
in  our  computer  library.  This  algorithm  uses  the  orthogonality  properties  ot 
the  Chekyshev  polynomials.  The  development  is  analogous  to  that  followed  in 
producing  the  Gauss -Legendre  and  Gauss-Laguerre  quadratures.  In  this  case, 
the  pertinent  integral  and  corresponding  Gauss-Chekyshev  quadrature  formula 
are  given  by 


1 

/ -J- 

-1  /77 


F(z)  dz  * Z w1  Ffzj) 
1*0 


The  Integration  Is  exact  If  F(z)  is  a polynomial  of  degree  2n  + 1 or  less. 

Note  the  singularity  at  z = t 1.  The  focal  position  for  a fan  of  input  rays 
Is  given  in  Table  1 . These  data  are  for  an  F/1.5  lens.  The  heights  and  focal 
positions  are  normalized  to  the  lens  radius  of  1. 
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Table  1 


Input  Height 

Focal  Position 

.00684 

.99962 

.04496 

.99977 

.11466 

.00002 

.16573 

.99999 

.21358 

.99996 

.25728 

.99996 

.29596 

1 .00001 

Previous  results  using  Simpson  rule  and  Gauss-Legendre  quadratures  were  at 
best  accurate  to  two  figures.  In  some  cases  the  results  were  off  by  factors 
of  2 or  3.  Presently  we  can  trace  rays  through  a geodesic  lens  with  at  least 
four  figure  accuracy.  This  is  felt  to  be  adequate  for  purposes  of  further 
tolerance  analysis,  system  design  and  evaluation  of  IOSA  concepts.  This 
program  does  not  require  the  equivalent  lens  transformation. 

It  is  also  possible  to  design  lenses  by  combining  this  new  capability 
with  existing  optimization  routines.  One  method  for  doing  this  is  to: 

(1)  Define  several  points  on  an  initial  profile 

(2)  Fit  a cubic  spline  curve  through  these  points.  (This  curve  has 
zero  slope  at  the  edge  and  center  of  lens.) 

(3)  Use  optimization  routines  to  determine  profile  by  varying 
the  depth  at  the  initial  points.  This  final  step  requires 
ray  tracing. 

(c)  Tolerance  Analysis 

Eq.  15  provides  a very  convenient  way  of  performing  tolerance  analysis 
of  a geodesic  lens  design.  Given  the  profile  and  the  variation  of  interest, 
the  slope  at  each  point  in  the  lens  can  be  determined.  In  a focal  system 
such  as  shown  in  Figure  11,  the  separation  between  a ray  and  the  optical 
axis  for  any  distance  from  the  lens  center,  f,  is  given  by 

d = R sin  0£  + (f  + R cos  6g)tan  (60  +6g)  (21) 
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of  the  geometrical  spot._siz&. — _The  focal  length  is  defined 
as  the  length  f where  4(<Mf))Z  is  a minimum.  The  spot  size 
is  the  full  extent  of  the  rays  at  this  point. 


*£  * 


The  distance  from  the  center  of  the  lens  to  the  ray  intersection  with  the 
optical  axis  is  given  by 


For  a lens  that  does  not  conform  to  the  optimum  profile,  F is  a function 
of  the  ray;  hence,  the  focal  position  for  the  tolerance  analysis  was 
chosen  as  the  position,  fm,  of  minimum 

0(f)  = A d.(f)  2 , (23) 

i 1 

where  d^(f)  is  the  distance  of  the  i^  ray  from  the  optical  axis  at  a distance 
f from  the  center  of  the  lens.  Figure  11  illustrates  these  parameters.  The 
minimum  is  found  by  numerically  searching  for  the  minimum  by  varying  the 
distance  f between  the  extremes  in  the  focal  position  of  the  individual  rays. 

Two  functional  variations  from  the  optimum  profile  have  been  analyzed. 

In  the  first,  the  lens  depth  was  varied  by 

Ay/y  = constant,  (24) 

where  Ay  is  the  variation,  and  y is  the  initial  depth.  Both  Ay  and  y are 
functions  of  the  radial  position,  p.  This  particular  variation  represents 
an  error  caused  by  tool  we ar  during  impact  grinding.  In  the  second  variation 
a percentage  of  the  upper  portion  of  the  lens  was  neglected  while  the  profile 
was  maintained.  Experimentally,  this  corresponds  to  either  not  grinding  deep 
enough  or  to  polishing  away  a portion  of  the  lens  rounding. 

The  profile  of  an  f/1.5  lens  was  used  in  the  calculations;  some  results 
are  given  in  Figs.  12throughI4.  Figure  12  shows  the  increase  in  focal 
length  versus  the  percentage  decrease  in  depth  for  tool  wear  and  for  surface 
removal.  For  each  unit  decrease  in  total  depth  for  tool  wear,  there  is 
approximately  a 13  unit  change  in  focal  length.  With  surface  removal,  the 


Focal  length  increase  vs  % of  depth  decrease 
for  both  tool  wear  and  surface  removal . 


A spot  diagram  at  the  focal  position  for  100  rays  and  a 2%  tool  wear 

depth  change  is  shown  in  Figure  14-  The  most  interesting  thing  about  this 

diagram  is  the  tight  core  and  the  few  scattered  rays  outside  of  this  core. 

These  scattered  rays  are  those  passing  through  approximately  the  outer  10% 

of  the  original  useful  aperture.  Another  point  of  interest  is  that  all 

the  spot  diagrams  provided  have  this  feature  of  a tighter  core  for  90% 

of  the  useful  aperture.  This  is  further  illustrated  in  Figure  14,  where 

the  geometrical  spot  sizes  for  full  aperture  and  for  90%  aperture  are 

plotted  versus  the  % depth  decrease. 

For  an  f/1.5  lens  in  LiNbO^,  the  full  lobe  width  of  a Gaussian  beam 

truncated  at  the  1/e  point  is  ^1.4  urn  at  a free  space  wavelength  of  0.8  (1 

-4 

ym.  For  a 1-cm-diameter  lens,  this  corresponds  to  3 x 10  units  in 
Figure  15.  Hence,  to  maintain  diffraction  limited  performance,  a depth 
change  of  1%  can  be  tolerated  if  90%  of  the  original  useful  aperture  is 
utilized.  For  even  worse  tool  wear,  a smaller  aperture  is  required  to 
maintain  diffraction  limited  performance. 


Fig.  13  Spot  diagram  for  a 2.0£  total  tool  wear  depth  decrease. 


SECTION  VI 


CONCLUSIONS  AND  SUGGESTIONS 


Under  this  program,  we  have  successfully  developed  the  ultrasonic 
impact  grinding  technique  for  the  fabrication  of  aspherical  geodesic  lens 
depressions.  The  lens  radius  is  4 mm  and  the  depth  is  1.374  mm.  15X  of 
the  lens  aperture  is  allocated  for  a curved  rounding  at  the  edge  of  the 
depression.  The  designed  focal  length  is  10.2  mm  for  an  aberration  free 
aperture  of  6.8  mm,  thus  the  f number  of  the  lens  is  1.5.  The  technique 
of  ultrasonic  impact  grinding  is  proved  to  be  precise,  fast  and  economic. 
Reproducible  results  can  be  easily  obtained.  The  tool  wear  problem  can 
be  neglected  if  more  than  two  tools  are  used  for  grinding  the  lens.  We 
have  demonstrated  that  depth  variation  within  7.5  pm  for  5 depressions  can 
be  obtained  using  two  tools.  One  of  the  5 depressions  has  the  exact 
designed  depth.  The  consistency  of  the  lens  depressions  can  be  further 
improved  if  a third  tool  is  used.  The  grinding  process  is  fast  using 
the  200  W Bullen  machine.  It  takes  only  45  sec.  to  grind  a 1.4  mm  deep, 

8 mm  diameter  depression  in  LiNbO-j  substrate.  The  machine  is  easy  to 
operate  once  the  optimum  grinding  parameters  are  set. 

Aberration  free  geodesic  lenses  have  been  fabricated  in  Ti:liNb03 
waveguides.  The  measured  focal  length  is  fairly  close  to  the  designed 
10.2  mm  even  though  the  lens  profile  is  slightly  off  the  ideal  profile.  The 
major  problem  encountered  is  the  polishing  of  the  lenses.  There  are  many 
fine  scratches  at  the  edge  of  the  depression,  causing  excessive  propagation 
loss  and  in-plane  scattering.  The  measured  3dB  focal  spot  size  is  1.6  pm 
as  compared  to  the  theoretical  value  of  0.32  pm.  The  discrepancy  is  believed 
to  be  the  sum  of  the  following  three  effects.  (1)  The  resolution  of  a 

perfect  lOx  microscope  objective  with  NA  * 0.3  is  0.93  pm.  (2)  The  use  of 
prism  coupler  to  extract  guided  light  out  of  waveguide  causes  aberration. 

The  rutile  prism  is  uniaxial  and  will  introduce  additional  aberration. 

(3)  The  actual  lens  profile  is  slightly  off  the  ideal  designed  profile. 

More  sophisticated  measurement  technique  has  to  be  developed  to  evaluate  small 
f number  lens  with  submicron  focal  spot  size.  The  polishing  process  also 
has  to  be  improved,  otherwise,  the  in-plane  scattering  due  to  scratches  will 
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bury  the  side  lobe  information  and  limit  the  useful  dynamic  range  of  the 
lens  system.  It  is  important  to  emphasize  that  the  polishing  process  should 
not  destroy  the  lens  profile. 

An  extensive  computer  error  analysis  is  completed  under  the  program. 

The  purpose  of  this  study  is  to  determine  the  effect  of  deviation  of  the 
actual  profile  from  the  ideal  profile  on  the  spot  size  and  focal  length. 

A ray  tracing  through  geodesic  lens  depression  is  developed  from  which  the 
geometric  spot  size  and  focal  length  are  calculated.  For  an  f/1.5  lens  in 
LiNbOj,  the  full  lobe  width  of  a Gaussian  beam  truncated  at  the  1/e  points 
is  about  1.4  ym  assuming  the  free  space  wavelength  is  0.8  ym.  For  a 1 cm 
diameter  lens,  a depth  change  of  1%  can  still  maintain  its  diffraction  limited 
performance  if  90%  of  the  original  useful  aperture  is  utilized.  On  the 
other  hand,  the  focal  length  is  very  sensitive  to  the  depth  variation.  We 
have  shown  that  for  each  unit  decrease  in  total  depth  for  tool  wear,  there 
is  approximately  a 13  unit  change  in  focal  length.  With  surface  removal, 
the  change  in  focal  length  is  about  3 units  per  unit  change  in  depth. 

Unless  one  can  control  the  depth  better  than  1.0  ym,  the  integration  of 
geodesic  lenses  in  integrated  optic  circuits  is  an  extremely  difficult  job. 

At  the  current  stage,  in  any  optic  system  using  geodesic  lenses,  the  precise 
location  (within  several  ym)  of  the  focal  plane  will  have  to  be  determined 
after  the  lens  is  fabricated. 
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